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ABSTRACT 

In the minimal supersymmetric Standard Model, the heaviest CP-even Higgs boson H 
and the CP-odd Higgs scalar A are predicted to be almost degenerate in mass at the tree 
level, for the wide kinematic range > 2Mz and tan/3 > 2. However, if large soft-CP- 
violating Yukawa interactions involving scalar quarks of the third generation are present in 
the theory, then the CP invariance of the Higgs potential can be maximally broken beyond 
the Born approximation, and the high degree of mass degeneracy between H and A may 
be lifted through a sizeable HA mixing. After taking the CP-odd tadpole renormalization 
of the A boson into account, we find that the small mass difference Mh — Ma, which is 
about 1% of the A-boson mass at the tree level, can be substantially enhanced to the 25% 
level at the one-loop order. We also find that the loop-induced mixing between the lightest 
CP-even Higgs boson h and the A boson may be of comparable size to Mh- We briefly 
discuss the main phenomenological implications of the predicted hA and HA mixings for 
the general Higgs-boson mass spectrum and for CP-violating observables at collider and 
lower energies. 
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It is known that the Higgs potential of the minimal supersymmetric Standard 
Model (MSSM) is invariant under the transformations of charge conjugation and parity 
(CP) at the tree level. The reason is that supersymmetry (SUSY) imposes an additional 
(holomorphic) symmetry on the Higgs sector of a general two-Higgs doublet model, which 
entails flavour conservation in tree-level neutral currents and absence of CP- violating scalar- 
pseudoscalar mixings in the Born approximation. Beyond the Born approximation, recent 
studies have shown that CP invariance of the Higgs potential may in principle be broken 
spontaneously through radiative corrections || if the CP-odd Higgs scalar A is sufficiently 
light mi]. However, this possibility has now been ruled out by experiment Hf§. 

Here, we shall study in more detail another interesting possibility of CP non- 
conservation within the context of the MSSM. In this case, CP violation in the Higgs 
sector is induced by loop effects due to the presence of additional interactions in the the- 
ory, which violate CP explicitly. For example, such CP-violating interactions may occur 
in the trilinear Yukawa couplings of the Higgs fields to scalar quarks. However, the phe- 
nomenological viability of such a scenario of explicit CP violation is often questioned in 
the literature. The standard reasoning for the latter goes as follows. Unless the scalar 
quarks have TeV masses, all new CP-violating phases in the MSSM not present in the SM 
fU||,|7j must be suppressed at least by a factor of order 10~ 3 , otherwise one-loop effects may 
exceed the current experimental limit on the neutron electric dipole moment (EDM). Ac- 
cording to these general arguments, CP violation in the Higgs sector of the MSSM, which 
arises at one loop, was estimated to be dismally small, so as to bear any phenomenological 
relevance ||. In this paper, we shall show that without any further assumptions on the 
model, this rough estimate based on the neutron EDM limit and the naive counting of 
loop suppression factors is fairly inaccurate and may be misleading in general. In partic- 
ular, we find that scalar-pseudoscalar mixings can still be large in the MSSM giving rise 
to observable CP-violating phenomena, if the Yukawa sector involving the scalar top and 
bottom quarks contains relatively large CP-violating couplings compatible with universal 
boundary conditions imposed by minimal supergravity models at the unification scale My. 

The MSSM introduces several new phases in the theory which are absent in the SM 
||. The mixing mass parameter /x involving the two Higgs chiral superfields in the super- 
potential, the soft-SUSY-breaking gaugino masses m\ (with A = g, W and B representing 
the gauginos of the SU(3) C , SU(2)^ and U(l)y gauge groups, respectively), the soft bilinear 
Higgs mixing mass m\ 2 (sometimes denoted as B/j, in the literature) and the soft trilinear 
Yukawa couplings Af are all complex numbers. If the universality condition at Mx is as- 
sumed, the gaugino masses mx are then related to each other and have the same phase, 
while the different trilinear couplings Af are all equal, i.e., Af = A. Not all phases of the 
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four complex parameters {//, mf 2 , m\, A} are independent of the fields ||. For instance, 
one can make m\ real by redefining the gaugino field A. Also, as we will see below, min- 
imization conditions on the Higgs potential lead to the constraint that the phase of m\ 2 
should be equal to the phase difference of the two Higgs doublets in the MSSM. As a re- 
sult, there are only two independent CP- violating phases in this constrained version of the 
MSSM. Usually, these are taken to be arg(//) and arg(A). Limits coming from the electron 
EDM may be avoided to a great extent by requiring that the phases of m\ and /i be aligned 
to Sm(m» = 0| . In this case, chargino and neutralino mass matrices are 

real and conserve CP. Furthermore, it has been argued |§ that bounds obtained from the 
neutron EDM leave arg(A) essentially unconstrained, even for slightly smaller than TeV 
soft scalar masses || . Notice that constraints on the scalar- top and bottom sector do not 
come directly from the neutron EDM but rather indirectly via the universality conditions 
at Mx- From the above discussion, it is clear that for soft scalar masses in the range 
0.5 < M < 1 TeV, a.rg(A) can safely be considered to be the only large CP-violating phase 
in the theory of order unity. As a consequence, significant scalar-pseudoscalar mixings in 
the MSSM are only expected to come from loop effects of scalar top and bottom quarks, 
whereas chargino and neutralino contributions being proportional to the vanishingly small 
CP-violating phase arg(/i) may be neglected. 

We start our discussion by considering the Higgs potential of the MSSM. Because of 
the holomorphic property of SUSY, one needs two Higgs doublets at least, denoted as $i 
and $ 2) with opposite hypercharges, K(<3> 2 ) = — = 1, in order to give masses to both 
up- and down- quark families and, at the same time, cancel the triangle anomalies induced 
by the fermionic SUSY partners of the Higgs field. After integrating over the Grassmann- 
valued coordinates in the SUSY action and including the soft-SUSY-breaking masses for 
the Higgs fields, one arrives at the Lagrangian describing the Higgs potential 

+ A 2 (<K 2 $ 2 ) 2 + A 3 (*I$i)(S^2) + A 4 ($ t 1 $ 2 )($ 2 <l> 1 ) , (1) 
where $i = — ir 2 $^ (r 2 is the Pauli matrix) and 

2 2 I 12 2 2 I |2 

fh = ~ m \ ~ m ' ^2 = ~ m 2 - W > 

Ai = A 2 = -Itf + f), A 3 = ~(g 2 -g' 2 ), A 4 = \g\ (2) 

The complex parameter m\ 2 m Eq. ([[]) as well as the real parameters m\ and m\ in Eq. 
(0) are soft Higgs-scalar masses. Furthermore, g and g' are the usual SU(2)^ and U(l)y 
gauge couplings, respectively. It is interesting to remark that the MSSM Higgs potential 
Cy contains fewer quartic couplings than that of the general two-Higgs doublet model; all 
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quartic couplings Aj in Cy are uniquely specified by SUSY. This makes the Higgs sector of 
the MSSM highly predictive. 

In order to determine the ground state of the MSSM Higgs potential, we consider the 
linear decompositions of the Higgs fields 

= ( * ) = # ( 02+ ) (3) 

\j ! (v l + H 1 + iA l ))' 2 [j. (v 2 + H 2 + iA 2 ) ) ' [6 > 

where V\ and v 2 are the moduli of the vacuum expectation values (VEV's) of the Higgs 
doublets and £ is their relative phase. The positive parameters V\ and v 2 and the phase 
£ are entirely fixed by the minimization conditions on Ly. This can be accomplished by 
requiring the vanishing of the following tadpole parameters: 



dC 



T Hl = (^) = vi (4 + Ke(m 2 2 e i€ ) tan/3 - ^ Mf cos 2(3 , (4) 



dC 



T H2 = (^) = v 2 nl + Ke(m 2 2 e i? ) cot (3 + £ M\ cos 2(3 , (5) 



dH 2 



2 
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= (^} = ^mKe^), (6) 

T A2 = (gf 2 ) = -^mK^l, (7) 

where tan (3 = v 2 /v± and M§ = (g 2 +g' 2 )v 2 /4 is the Z-boson mass squared with v 2 = v\+v 2 . 
Variations of Cy with respect to <pf and <p 2 vanish identically, reflecting the fact that a 
physical ground state should conserve charge [f|. If we now perform the orthogonal rotation 




cos (3 — sin (3 
sin (3 cos (3 




(8) 



the Higgs potential shows up a flat direction with respect to the G° field, i.e., (dCv /dG°) = 
0. Then, the G° field becomes the true would-be Goldstone boson eaten by the longitudinal 
component of the Z boson. In this weak basis, the tree-level mass matrix of the CP-odd 
scalars becomes diag(0, M|), where M\ = 9?e(m 2 2 e^) / '{spcp) is the tree-level A-boson 
mass squared. Furthermore, the orthogonal rotation @ leads to a non-trivial CP-odd 
tadpole parameter given by 

Ta = { ~M ] = - v ^ m ( m i2^)- ( 9 ) 



As has been shown explicitly in fLOl , the tadpole renormalization of the A boson is very 
crucial to render all HiG°, H 2 G°, H\A and H 2 A mixings ultra-violet (UV) finite (see also 
Fig. 1). It is now important to notice that the phase difference £ between the two Higgs 
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VEV's is no more arbitrary but completely specified by virtue of Eq. (0). At the tree level, 
one has T4 = and m\ 2 e 1 ^ is a real number. Beyond the Born approximation however, 
m 2 2 e acquires a small imaginary part. In fact, the no n- vanishing tadpole graph of the 
A-boson should be compensated by the tadpole counter-term (CT) Ta, such that the true 
ground state of the effective Higgs potential does not get shifted. Without any loss of 
generality, we can therefore redefine $2 as e~* arg( - m i 2 - ) $2 ) thereby resulting in a weak basis 
in which £ = at zeroth order. Then, m\ 2 ls rea l a ^ the tree level but effectively receives 
an imaginary part at higher orders which is determined by the tadpole renormalization 
condition on Ta- With this simplification, one can avoid unnecessary ^-dependent phases 
in the MSSM. 

After including all tadpole contributions, the Lagrangian relevant for the Higgs-boson 
mass matrix can be cast into the general form 



C 



H 



Hi, H 2 , G°, A 



Mis 



Ml 



Employing the usual short-hand notations s 3 
in Eq. ([10|) are given by 



sin x and cos x 



(10) 



/ Hi \ 
H 2 
G° 
V A J 
the 2x2 sub-matrices 



M% 



Ml 



M 2 SP 



cj|Mi + s%M\ - T H Jv x 
-s p c p (M 2 z + Mi) 



spT Hl — cpT H2 
v 



S/3 Cp 
-C/3 S/3 




Mi 



-spepiMi + Mi) 
slM% + clM\ - T H Jv 2 



spT Hl 



- cpT H2 



S/3 tan (3 Th 1 + cp cot f3 Th 2 
v 



/ 



(11) 



(12) 



(13) 



and M 2 PS = (M SP ) T . In the Born approximation, the CP-even mass eigenstates h and H 
are obtained by diagonalizing M 2 S through the orthogonal transformation 




cos a 
sin a 




(14) 



with 



tan(2a) 



The tree-level mass eigenvalues of M\ are then given by 

1 
2 



Ml 



h(H) 



M 2 z + M 2 A 



(M 



AM\M 2 Z cos 2 2/3 



(15) 



(16) 
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Notice that for tan/3 > 2, one has cos 2 2(3 ss 1, M h ps M 2 and M# ps Ma, that is, if and 
A are nearly degenerate in the MSSM. However, radiative corrections due to the large top 
Yukawa coupling affect the /i-boson mass in a significant manner. In fact, h can be as heavy 
as 130 GeV and is heavier than the Z boson for the largest bulk of the parameter space 
|TT| . On the other hand, the high degree of mass degeneracy between H and A generally 
persists even beyond the tree level in the CP-invariant limit of the theory, especially when 
Ma > 2Mz and tan (3 > 2. In this kinematic range, tan (3 ~ tan a and quantum effects seem 
to affect equally Mh and Ma, such that the small mass splitting of H and A remains still 
valid [p!2l . As we will see however, large CP- violating Yukawa couplings of scalar top and 
bottom quarks can give rise to sizable HA mixings at the one loop level and to a substantial 
enhancement of the mass difference Mh — M&. In the calculation of scalar-pseudoscalar 
self-energy transitions, one has to include the relevant CP-violating mass CT's given by 
Ai'sp in order to arrive at UV-safe analytic results. 

We shall now discuss the interactions of the neutral Higgs fields with the scalar top 
and bottom quarks in the presence of CP- violating terms [13]. The scalar quarks of the 
third generation are represented by 



Q = r h = t R , D* = b R , (17) 




with Y(Q) = 1/3, Y(U) = —4/3, and Y(D) = 2/3. If we denote the neutral component of 
the Higgs doublets $i and $ 2 by cf)\ and there are then three contributions of <p\ and 
02 to scalar-quark masses and their respective couplings. The relevant Lagrangian, Cy, 
receives contributions from the soft-SUSY-breaking sector, and from the so-called F and 
D terms (or components) of the chiral and vector superfields, respectively. Specifically, we 
have 

£y = ^soft + + £d > ( 18 ) 

where 

-C° so{t = M 2 Q (tlt L + blh) + M?t* R t R + Mfb R b R + (fiA b ^b R b L + f 2 A t <p° 2 t R t L 
+ H.c.), 

-£% = AVM(^l + b* R b R ) + fUl*&{ili L + t R i R ) 

+ (V2 4>°x%t R + nh 4>° 2 b* L b R + h.c. ) , 

- C% = -2- - 02>°) [(1 - 2e,4) ilh + 2e t s 2 w t*J R - (1 + 2e b s 2 J b* L b L 

+ 2e b s 2 J R b R ], (19) 
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with s w = sm9 w being the weak mixing angle, e t = 2/3, e& = —1/3, fi = v2m&A'i and 
f2 = y/2m t /v2- Furthermore, Mq, M t and Mb are soft-scalar quark masses and are usually 
considered to be all equal to M at the unification scale Mx- Notice that all operators of 
dimension four in Cy and C v satisfy an additional global U(1)q symmetry with Q charges: 
Q($i) = 2, Q($ 2 ) = 1, Q(Q) = 0, Q(U) = -2 and Q(D) = -1; the bilinear operator $1$ 2 
and the A-dependent trilinear terms break it only softly. In fact, if the Higgs- mixing term 
(jl in the superpotential is absent, the U(1)q symmetry can appropriately be extended to 
the whole MSSM Lagrangian. As has extensively been discussed in and is also valid 



for the case of the MSSM, the simultaneous soft-breaking of the symmetries U(1)q and CP 
is sufficient to assure the renormalizability of Higgs scalar-pseudoscalar transitions to all 
orders in perturbation theory. 

It is straightforward to obtain the scalar top and bottom mass matrices from C v 
which may conveniently be expressed as 



jtf = [ Ml + m\ + cos2/?Ml (T| - e q s 2 w ) m q (^R q + A*) e* - , 

' m q (fJ l *R q + A q ) e~ iSq M 2 + m\ + cos2/3Mf e g s? 



where q = t,b, T\ = 1/2, T\ = —1/2, R t = cot ft, Rb = tan ft. The phase S q is determined 
from the requirement that the scalar-quark mass matrix becomes positive definite by a 
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judicious re-definition of the right-handed scalar quark fields, i.e., q^ — > e q qn with 5, 
axg(fi*R q + A q ). In this weak basis, we can diagonalize A4 2 through the orthogonal rotation 

(21) 

COSUq I \ q 2 I 

where the rotation angle 9 q may be obtained by 




cos 6q sin 6 q 
— sin 6q cos 9 q 




2m a \uR + A*\ , , 

tan (20 p ) = - * ^—^ q - . 22 

V " ; M 2 -M$ + cos2ftM 2 (T«-2e q s^ 



w J 



The masses for two scalar-quark mass eigenstates qi and g 2 are then given by 
1 



M m.) = 2 \ M Q + M ,+ 2m q + T ^os2ftM' z 



+ 



M 2 - Ml + cos2ftM 2 (Tl - 2e q s 2 w )] 2 + 4m 2 \iiR q + A* q \ 2 \. (23) 



It is easy to see that the scalar quarks of the first two families are almost degenerate for 
large universal soft-scalar quark masses, e.g., bigger than 0.5 TeV. However, this mass 
pattern is not in general true for the third family owing to the non-negligible top and 
bottom masses. 
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Defining as ip° k = a/2 (0° — (0°)) = H k + iA k with k = 1, 2, we are then able to write 
the trilinear couplings of the Higgs fields to the scalar quarks in the generic form 



C 



hit 



q=t,b 



L,q , L,<3 \ ~* ~ 



R.q , -R,g \ ~* ~ 



where the different coupling parameters are given by 



(24) 



L,t 

9i 



R,t 

9i 



L,b 

9x 



R,l 

9i 



V 



- e t s, 



I 2 



#2 



mi 



+ 



v ^2 



C/3 M 



(1 + 



92 



mr 



_ sjMl A _ , x 

7, V9 



L,6 
#2 



v ^2 



(25) 



b _ m & A*„-i8 b 



K = —Ale 



m b 



h 2 = \xe 



-I0 b 



From Eqs. fl24]) and (p5|), it is then not difficult to compute the tree-level couplings Hiq*qj, 
H 2 q*qj, Aq*qj and G°q*qj = 1,2). These couplings may respectively be given by the 
following matrices: 



^H k q*q 




^Aq*q 



s 2q Ueh q k - c 2 q g k ' q - s 2 q g k ' q -c 2q Steh\ - s q c q (g k ' q - g k ' q ) + 2iT^mh q k 
-c 2q ^eh q - s q c q (g k ' q - g k ' q ) - 2iT^mh q k -s 2q Reh q k - s 2 q g k ' q - c 2 q g k ' q 



i 



( 



T%(meh\ - 



—is 2q ^smh\ 
cps 2q AM 2 q 



-T*(23te/i? 



c p s 2q AM 2 



+ ic 2q ^smh\ 



+ ic 2q ^smh\ 



is 



2q ^smh\ 



J 



J- r 



T. 



q is 2q AM 2 ~ 



1 
-1 



(26) 



where k = 1,2 and AM? = M| - M| r Note that SmftJ = - tan/5 3m/i^ and Ueh\ = 
— tan/? 3?e/i2 + ( s q c qAM~ /v) in the weak basis we are working. 

As has been discussed above, the tadpole parameter T4 is determined by the A- 
boson tadpole graph, T A (0), shown in Fig. 1(c). Using the renormalization condition 
T A + T A (0) = 0, we easily find that 

d n k 



- E nz 

q=t,b 
1 



(2n) n i 

S2q 



Tr 



iA q (k)iT, 



Aq'q 




167T 2 



E K ^ <Zmh\ AM- So(0,M?,M|) 



q=t,b 



80 



(27) 



S 



In Eq. = 3 is the colour factor for quarks, A«(Jfe) = diag[(fc 2 -M?)" 1 , (P-M?)- 1 ] 

is the scalar-quark propagator matrix and B (s, mf, mf) is the known Veltman-Passarino 
function defined as 



C 



uv 



ln(mim 2 ) + 2 + - (ml — m?) In f — - ) 

s 1 v ra2 y 



+ A 1 / 2 ( S ,m 2 ,m 2 ; 



mi + ro. 



2m\m2 



(2? 



where A(x, ?/, 2) = (x — y — z) 2 — Ayz and C 



uv 



1/e — 7s + ln(47r/i 2 ) is an UV infinite 
constant. For s — 0, the one- loop function in Eq. (p8| ) takes on the simple form 

B o (0, m 2 ,m 2 ) = C uv - H mi m 2 ) + 1 + m \ + ^§ In ( — ) . (29) 

A straightforward calculation of the individual contributions to H^A self-energies (with 
% = 1,2), shown in Figs. 1(a) and (b), yields 

d n k 



q=t,b 



Tr 



iA q (k) iT^ r<} iA q (k + p) iTQ 



Aq'q 



167T 



(27r) n i 

l — Y, N? Zmhl { ^ AM- B^MlM 



q=t,b 



S/3 



£2g 

S /3 



+ ^f' 9 ) (*o(s, Ml Ml) - B (s, Ml Ml 



+ (s 2 q 9-' q + c 2 q g^) (B (s, Ml Ml - B (s, M 2 , M? )) ] } 
-(M 2 P h = ^ 



-L 2 K Zmhl ™* AM 2 i?o(0,M?,M 2 

167T 2 ^ ^ 



52/ ' 



(30) 



(31) 



where s = p 2 , r\ = cos/3 and r 2 = sin/5. It is easy to verify that the HiA self-energies 
are UV finite, only after the CP-odd tadpole contribution given in Eq. ([H]) is taken into 
account |TD] . The CP- violating HA and HA self-energy transitions are then obtained by 

(32) 



Ii hA {s) 
U HA (s 



s) + sin a Tl H2A (s) 



cos a Ii HlJ 

— sin a + cos a 



For the kinematic range Ma > 2Mz and tan/5 > 2 (tana ~ tan/5), we have H hA (s) ~ 
II i?2j4 (s) and Il HA (s) ~ — II i?lj4 (s). Furthermore, the rotation angles t and are almost 
equal to 7r/4, i.e., s 2g ~ 1, since the off-diagonal elements of the scalar-quark mass matrix 
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M. 2 q are much bigger than the difference of its diagonal entries. For relatively low tan/3 
values, the biggest contribution occurs when the mass difference between t\ and t<i is rather 
large. Specifically, for 2 < tan/? < 5, we have 

II feA (0) w 4^1 lUeh'J^^- \nx t ~ l) - ^MnzJ, (33) 
Sir 2 1 y xf — I ' v J 

n HA (0) « --^^^(4^111^-1), (34) 
oir z ^ xf — 1 ' 

with x 4 = MijMi x > 1. In Eqs. (|33|) and fl3~4[) , the CP- violating quantities ^smh\^R.eh[ and 
Qmhffieh!^ satisfy the inequalities 

1 777 ^ 1 777 ^ 

ISm/i^e/i'l < — ^ |/x| 2 , ISmtffte/i'l < - — L. I^J . (35) 

2 irsi 2 ?rsi 

For larger tan/? values, scalar-pseudoscalar mixings receive significant contributions from 
scalar-bottom quarks as well, leading to a more involved kinematic dependence. 

We shall now examine numerically the dependence of HA and hA mixings on the 
various kinematic parameters. Motivated by minimal supergravity models, we consider 
the soft scalar-quark masses and trilinear couplings to be universal: Mq = M t = M& = 
Mq = 0.5 TeV and A t = Af, = A with arg(A) = 90° |jl4] |. The soft parameters m\ and 



?t?| in the Higgs potential (c.f. Eq. (0)) are also equal to Mq at M x . However, their 
renormalization-group-equation (RGE) runnings from Mx down to Mz are very sensitive 
to further model-dependent details of the physics near the unification scale |IJ|. Therefore, 



we assume that m\ and m\ are free parameters not larger though than few TeV, being 
subject into the tadpole constraints in Eqs. @) and (||). For our numerical analysis, we 
consider the parameters {tan/3, /i, M , A} as independent and take Ma = M . We expect 
that the above considerations will be in good agreement with a more elaborate treatment. 
Refinements due to the RGE running of the MSSM parameters are beyond the scope of 
the present paper and may be studied elsewhere. 

In Fig. 2, we display the dependence of the HA and hA self-energies on the param- 
eter /i at vanishing momentum transfer, s = 0, for three selected values of \A\. With \i 
increasing, both the CP-violating quantities in Eq. ( P5[ ) and the scalar-top mass ratio x t 
increase and so the HA and hA mixings. Indeed, we observe that the mass difference 
Mu — Ma ~ U HA (0) /M A can be as large as 25% x M A , for /i ~ 2 TeV (M A = 0.5 TeV) and 
Mj i5 M^ > Mz- This should be contrasted with the tree-level prediction of 1.07% x Ma 
for the same value of tan/3 = 2. As has been discussed in |lT0|Jl~5[] , an HA mass splitting 
of order 10% of the A-boson mass is sufficient to give rise to significant contributions to 



the neutron EDM through the two-loop Barr-Zee mechanism |T6| close to the observable 
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limit | I7[| . Furthermore, we find that U hA (0) may reach the (100 GeV) 2 level in the same 
kinematic region and so be comparable in size with at one loop. Nevertheless, the net 
effect of hA mixing on the /i-boson mass will be a modest reduction of the Mh one-loop 
value by -(n feA ) 2 /(2M|M|), ranging from -2% to -10% (M h » 100 GeV) in the A-boson 
mass interval 500 - 200 GeV. 

In Fig. 3, we plot the dependence of the scalar-pseudoscalar mixings as a function of 
energy (E cms = y/s). We find that — 3?en hj4 (s) (solid line) can increase by one order of 
magnitude for s ~ M\ in comparison to the s = value. The cause for this strong energy 
dependence appears to be the subtle cancellations occurring at low energies between terms 
proportional to dteh^ and m 2 /v in Eq. (|33| ) which get less important at higher energies. For 
similar reasons, $teIl HA (s) (solid line) shows up a milder energy dependence because of the 
absence of such destructive terms (see also Eq. fl3~ID). As can be seen from Fig. 3 too, the 
absorptive parts QmIl HA (s) and QmIl hA (s) (dashed lines) become significant above the i\ti 
threshold. In general, we have that both |3?eII- ffyl (M^)| and \QmU HA (M\)\ can formally 
be of comparable order. Thus, the necessary conditions for resonant CP violation through 
HA mixing |18| at future high-energy machines, such as the LHC, / u + collider, etc., may 



comfortably be satisfied ||15|| . In particular, at a \i \i collider with an integrated luminosity 



of 50 fb 1 , one may be sensitive to HA mixings up to order 10 3 TeV 2 , corresponding to 
mass differences M H - M A ~ 1% x M A . This is in agreement with an earlier observation 



made in |L8f that CP asymmetries of order unity due to resonant Higgs scalar-pseudoscalar 
transitions can naturally occur even within the MSSM |[L9| . Finally, the dependence of HA 
and hA mixings on tan/5 is presented in Fig. 4. We see that U HA (0) decreases in general 
for large values of tan/3, while — U hA (0) exhibits a significant enhancement for tan/? ~ 50 
close to the perturbative bound of the theory. The reason for this enhancement is due to 
the sizable left-right mixing of the scalar bottom quarks in the large tan (3 domain. 

The fact that scalar-pseudoscalar mixings may be large within the MSSM should 
not be very surprising and may be understood in simple terms as follows. All quartic 
couplings of the MSSM Higgs potential as well as Mh are suppressed by gauge coupling 
constants as a result of SUSY. On the other hand, despite the typical loop suppression 
factor (47r)~ 2 , Higgs scalar-pseudoscalar mixings are enhanced by the large top Yukawa 
couplings proportional to rut/i/v 2 and m t A t /v 2 . In this case it is therefore evident that 
naive dimensional analysis of counting loop suppression factors can lead to a dramatic 
underestimation of the actual size of the quantum effects. As has been shown by the present 
paper, a large Higgs scalar-pseudoscalar mixing is indeed possible within the MSSM and 
can lead to observable CP-violating phenomena in near future experiments. 
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Figure captions 



Fig. 1: Feynman graphs contributing to the Hi A and H 2 A mixings: (a) One-loop 
self-energy graph, (b) CP-odd tadpole renormalization, (c) Tadpole graph of 
the A boson. 

Fig. 2: HA and hA self-energies at s = as a function of the parameter /i. 

Fig. 3: HA and hA self-energies versus energy E cms = yfs. Solid lines indicate 
dispersive parts of the scalar-pseudoscalar self-energies, while dashed lines 
correspond to their absorptive parts 

Fig. 4: HA and hA self-energies at s = as a function of tan (3. 
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